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In this paper, we report the phenomena of global and partial phase synchronizations in linear 
arrays of unidirectionally coupled piecewise linear time-delay systems. In particular, in a linear 
array with open end boundary conditions, global phase synchronization (GPS) is achieved by a 
sequential synchronization of local oscillators in the array as a function of the coupling strength 
(a second order transition). Several phase synchronized clusters are also formed during the 
transition to GPS at intermediate values of the coupling strength, as a prelude to full scale 
synchronization. On the other hand, in a linear array with closed end boundary conditions 
(ring topology), partial phase synchronization (PPS) is achieved by forming different groups of 
phase synchronized clusters above some threshold value of the coupling strength (a first order 
transition) where they continue to be in a stable PPS state. We confirm the occurrence of both 
global and partial phase synchronizations in two different piecewise linear time-delay systems 
using various qualitative and quantitative measures in three different framework, namely, using 
explicit phase, recurrence quantification analysis and the framework of localized sets. 

Keywords: Global phase synchronization; partial phase synchronization; piecewise linear time- 
delay systems. 



1. Introduction 

Chaotic phase synchronization (CPS) associated with a locking of the phases of coupled chaotic sys- 
tems corresponds to the case where th e amplitudes ar e still uncorrelated. CPS has been investigated 



in ensembles of globally coupled arrays [Pikovskv e^ al. j_ 2001: Boccaletti et al. , 200d: Ivanchenko et al. 



2004l:lTakamatsu et alJ. l2000l: iPikovskv e^ aLl.ll996l:lKiss et ali l200ll : IZhou et al.\ . l2002l:IOsipov et al 



Zhan et all l2000l : iKozyreff et all l2000l : lOtsuka et all l2006f | , networks of oscillators [Boccaletti et al. 
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Arenas et al. 

2nnnl: 



200^ 



Otsuka et al. 



Bartsch et all 120071] . ecology 



Batista et ail . boOTJ : iRen fc Zhaol . |2007I:IYu et aUl2009ll. lase r systems iKozvreff et al. 



20061 ] , c ardiorespirato r y systems [Scha 



matology [Rybski et all l2006l : 



Yamasaki et al 



Blasius et al . 19991: Sismondd. 



2009[ 



Maraun et al 



1990l : 

?t all 



'er et all Il998l: IStefanovska et al 



Amritkar fc Rangaraianl . l2006l |. cli- 



2000 



20051 ] ■ etc. CPS is well studied and 



understood in low dimensional systems; however, there exist a very little indepth studies in higher dimen- 
sional systems such as time-delay systems, which are essentially infinite-dimensional and exhibit highly 
non-phase-coherent hyperchaotic attractors with complex topological structures. So estimating the phase 
explicitly to identify phase synchronization in such systems is quite difficult. Recently, the occ urrence of 
phase synchronization in time-delay systems has been reported Senthilkumar et ad . l2006l . l2008l ] . However, 
these investigations are carried out so far only in a system of two coupled time-delay systems. Very recently, 
the occurrence of global phase synchronization in a linear array of Mackey-Glass time-delay systems in 
the chaotic regime has been reported in Suresh et al, [2010] for open end boundary conditions. It has 
been shown that global phase synchronization occurs via a sequential type synchronization as the coupling 
strength increases. To verify the generic nature of the results of Suresh et al, [2010], we consider in this 
paper two different piecewise linear time-delay systems with complex topological structures and investi- 
gate the occurrence of global phase synchronization along with the underlying mechanism for open end 
boundary conditions. In addition we consider the case of closed end boundary conditions (ring topology) 
and identify the occurrence of partial phase synchronization via cluster formation. 

Specifically, in the first part of this paper, we investigate the generic nature of the phenomenon of 
global phase synchronization in an linear array, with free ends, of unidirectionally coupled (i) piecewise 
linear and (ii) threshold nonlinear time-delay systems in hyperchaotic regimes. In the second part, we 
report the phenomena of partial phase synchronizations in the array of above two systems but with closed 
end boundary conditions. At first, we use the nonlinear transformation introduced in Senthilkumar et al, 
[2006, 2008] to estimate the phase of all the systems in the array and to identify the existence of phase 
synchronization. Further, we confirm the existence of global phase synchronization (GPS) and partial phase 
synchronization (PPS) from the original non-phase-coh erent hyperchaotic attractors using two inde pendent 



Romano et all l2005l : iMarwan et all 120071 ] and the 



approaches, namely recu rrence quantification analysis ^ ^ 

concept of localized sets [Pereira et ad . 120071 ]. In addition, we point out that the onset of GPS in the linear 
array with open end boundary conditions takes place in the form of a sequential synchronization (second 
order transition) : For lower values of coupling strength the phases of nearby systems get already entrained 
with the drive system in contrast to the far away systems, while the other non-synchronized systems display 
clusters of phase synchronized states among themselves before they become synchronized with the large 
cluster in the sequence as the strength of coupling increases to form the GPS. 

On the other hand, if we consider an array with closed end boundary conditions (ring topology), partial 
phase synchronization occurs with the formation of different groups of phase synchronized clusters. The 
oscillators in the array self-organize to form groups of phase synchronized clusters where each cluster is 
in perfect phase synchrony as a function of the coupling strength. Such a clustering is considered to be 



m periect pnase syncnrony as a lunction or tne coupling strengtn. bucn a, clustering is considereu to be 
particularly significant in bi ological systems Kanekd . Il990l : IShermanl . Il994l : IStrogatz Sz Stewart! . 119931 ] , in 



electro-chemical oscillators [Kiss et all 1200111. etc. Re c ently, cluster synchronization in an array of three 



chaotic lasers without delay was reported [Terry et all 119991 ] as well. In our case, as the coupling strength 



increases, every individual system starts to drive the nearest system and the systems with small differences 
in their phases organize themselves to form small groups of clusters leaving the other oscillators with large 
phase differences to evolve independently. Further increase in the coupling strength results in an increase 
in the sizes of the clusters due to the locking of the phases of the nearby oscillators of the clusters for 
appropriate values of the coupling strength and finally ending up with large groups of clusters resulting in 
PPS (first order transition). 

The paper is organized as follows: In Sees. [3] and [H we will describe briefly the coupling configuration 
and the occurrence of global phase synchronization in an array of two different piecewise linear time-delay 
systems with open end boundary conditions. In Sees. [5] and [6] we consider the linear array of two piecewise 
linear time-delay systems with closed end boundary conditions (ring topology) and discuss the occurrence 
of partial phase synchronization and finally we summarize our results in Sec. [3 
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Fig. 1. Attractors of the piecewise linear system, (a) non-phase-coherent hyperchaotic attractor of the drive system (|4^). (b) 
Transformed attractor of the drive system, (c-f) Transformed attractors of some of the uniformly selected response systems 
{i = 10, 20, 30, 50) in the new state space {xi{t + r), Zi{t + r)). 

2. Estimates of phase 

Identifying Chaotic phase synchronization is a nontrivial problem in time-delay systems exhibiting com- 
phcated hyperchaotic attractors. For this purpose one requires appropriate measures. In this section, we 
describe briefly the various measures which we used from three different frameworks in this paper to es- 
timate the phase of time-delay systems thereby facilitating the characterization of phase synchronization 
between the coupled time-delay systems. 



2.1. Transformation of the original non-phase- coherent attractor 

Time-delay systems usually exhibit a highly complicated non-phase-coherent chaotic/hyperchaotic attrac- 
tor with their flows having more than a single center of rotation. Such an attract or does not yield a 



monotonica lly increasing phase on estimating it using the conventional techniques [Senthilkumar et al. 



20061 . I2OO8I ] ■ For this purpose, we have introduced a nonlinear transformation to recast the original non- 
phase-coherent attractor into a smeared limit cycle- like at tractor with a single center of rotation. Our 
transformation is obtained by defining a new state variable Senthilkumar et all l2006l . l2008l |. 



z{t + t) = x{t)x{t + f)/x{t + r), 



(1) 



where f is the optimal value of time-delay to be chosen in order to avoid any additional center of ro- 
tation. Then, the projected trajectory in the transformed state space {x{t -\- T),z{t + r)) will resemble 
that of a smeared limit cycle-like attractor. Such an approach facilitates the estimation of phase from the 
conventional techniques. 
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2.2. Recurrence quantification analysis 

Several measures of complexity quantifying the small scale structures in recurrence plots hav e been proposed 



and t he corresponding description is known as recurrence quantification analysis (RQA) [Marwan et al. 



20071 ] . In addition to the other advantages of the RQA such as its application to short experimental data, 
nonstationary and noisy data, we find that it can also be applied directly to highly complicated non- 
phase-coherent attractors of the time-delay systems in characterizing the synchronization transitions, in 
particular phase synchronization (PS) transitions. Among the available recurrence quantification measures, 
we use the Correlation of Probability of Recurrence (CPR) and the generalized autocorrelation function 
P{t), which are estimated from the original non-phase-coherent hyperchaotic attractors, to confirm the 
existence of GPS in the array of coupled time-delay systems, both qualitatively and quantitatively. 

A criterion to quantify phase synchronization between two systems is the Correlation of Probability of 
Recurrence (CPR) defined as 

CPR={Pi{t)P2{t))/aia2, (2) 

where P{t) is the recurrence-based generalized autocorrelation function defined as 

N~t 

^(*) = w~t ^ ®^' " "^^ " ^^+*"^' ^^^ 

where G is the Heaviside function, Xj is the i*^ data point of the system X, e is a predefined threshold, 
||.|| is the Euclidean norm, and N is the number of data points, Pi^2 means that the mean value has 
been subtracted and a"i^2 are the standard deviations of Pi{t) and P2{t), respectively. Looking at the 
coi ncidence of the positions of the max i ma of P{t) of the coupled systems, one can qualitatively identify 
PS Romano et al.l . l2005l : iMarwan et a?.l . 12007 ] . If both systems are in CPS, the probability of recurrence is 



maximal at the same time t and CPR ~ 1.0. If they are not in CPS, the maxima do not occur simultaneously 
and hence one can expect a drift in both probabilities of recurrences resulting in low values of CPR. 

2.3. The concept of localized sets 

Another measure that we have employed is the framework of localized sets Pereira et ali . 120071 ]. The 



concept of localized sets opens up the possibility of an easy and an efficient way to detect CPS especially 
in complicated non-phase-coherent attractors of time-delay systems. The main ingredient of this technique 
is to define an event in one of the systems and then observe the other during the event, which defines a set 
D. Depending upon the property of this set D, one can state whether PS exists or not. The set is spread 
over the entire attr actor for asynchronous systems, whereas it is localized on the attractor if the coupled 
systems are mutually phase-locked or in phase synchronous state. 

We use the above three different frameworks to confirm the existence of GPS in an array of piecewise 
linear time-delay systems in the following. 

3. GPS in a Linear Array of Time-delay Systems with piecewise linearity 

Recently, we have reported the dynamical organization of an array of Mackey-Glass time-del ay systems in 



chaot ic regime to form global phase synchronization (GPS) via sequential synchronization [Suresh et al. 



20101 ]. as ment ioned in the i ntrod uction. In this section, we intend to examine the generic nature of the 



results in Ref. Suresh et a/.L 120101 ] . by investigating the emergence of GPS in an array of a particular type 



of piecewise linear time-delay systems from the local sequential synchronization as the coupling strength 
is increased from zero to that of the hyperchaotic regimes. It is worth to emphasize that the hyperchaotic 
attractors of the piecewise linear time-delay systems are characterized by much more complex topo logical 



properties than the chaotic attractor of the Mackey-Glass time-delay system studied in [Suresh et aLl . l201Q] 



and so is of considerable significance in estimating the onset of GPS through sequential synchronization. 
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Fig. 2. Plot of the first ten maximal Lyapunov exponents An 
parameter values /3 = 1.0, ai — 1.2, r £ (3, 25). 



of the piecewise linear time-delay system (|4^) for the 



3.1. Linear array of piecewise linear time-delay systems 

Now, we consider an array of unidirectionally coupled piecewise linear time-delay systems with open end 
boundary conditions and with para meter mismatches. Linear stab i lity a r id bifurcation analysis of this 
syste m has been well studied in Ref. [Lakshmanan &: Senthilkuman . l20ld : ISenthilkumar Sz Lakshmanan . 
20051]. Many kind of sync hronizations and their transitions ha.ve been reported in the coupled piecewis e 
linear time-delay systems [Senthilkumar Sz Lakshmananl . l2005l . 120071 : iLakshmanan &: Senthilkuman . l2010l | . 
The array of unidirectionally coupled piecewise linear first-order delay differential equations is represented 



as 



xi{t) = -/3xi(t) + aif{xi{t - r)), 

Xi{t) = -f3xi{t) + aif{xi{t - t)) + e(xj_i(t) - Xi{t)), i 



2,3,--- ,A^, 



(4a) 
(4b) 



where a, /? are the system parameters, r is the time-delay and e is the coupling strength. We use the open 
end boundary conditions xq = xi and x^v = x^-i- The nonlinear function f{x) is chosen to be a piecewise 
linear function defined as 



/(^) 



0, x< -4/3 

-1.5a; -2, -4/3 <x< -0.8 

X, -0.8 < X < 0.8 

-1.5X-F2, -0.8 <x <4/3 

0, x> 4/3. 



(5) 



The system parameters are fixed as follows: /3 = 1.0, ai = 1.2, r = 15.0. The values of the nonlinear 
parameter Oj of the response systems in the array are chosen randomly in the range ai € (1-19, 1.23), so 
that all the subsystems are effectively nonidentical due to the parameter mismatch. The original non-phase- 
coherent hyperchaotic attractor of the drive in the {x{t),x{t — r)) state space is illustrated in Fig. \V[a). 
The corresponding transformed attractor, effected using the transformation ([T]), in the new state space 
{x{t + r), z{t + r)) is depicted in Fig. [T] (b), which now looks like a smeared limit cycle-like attractor with 
a single center of rotation. The optimal value of the f in ([T]) for the above piecewise linear time-delay 
systems is found to be 1.6 Senthilkumar et all 120061 ]. The same transformation has also been performed 
for the n qn-phase-coherent c haotic attractor of the Mackey-Glass time-delay system with an appropriate 
f in Ref. Suresh et al\ . \201u!\ . The hyperchaotic nature of the attractor in Figs. [1] for the above parameter 
values is confirmed from three positive Lyapunov exponents at r = 15.0 in the spectrum of ten maximal 
Lyapunov exponents of (|4^) as a function of the time-delay parameter r G (3, 25) in Fig. [2l 
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Fig. 3. (a-d) Plots of the phase differences {Acpii — (j>i 
for different e. 



5000 10000 15000 20000 
t 



of some of the selected systems (i = 10, 20, 30, 50) in the array 



4000 



2000 




e = 0.0 
£ = 0.56 

8 = 0.77 
£ = 0.79 





5000 10000 15000 20000 

t 

Fig. 4. The evolution of the average phase difference [ri{t)] for different values of the coupling strength e. 



3.2. Identification of GPS from the transformed attractor of the piecewise 
linear time- delay system 

To confirm the existence of GPS in the array of piecewise linear time-delay systems, we have fixed size of 
the array as A^ = 50. The transformed smeared limit cycle-like attractors with a single center of rotation of 
some of the uniformly selected subsystems (z = 10, 20, 30, 50 ) in the linear a r ray are shown in Figs.|T](c-f) . 
We have used the conventional Poincare section technique [Pikovsky et all I2OOII : iBoccaletti et all l2002l | 
to estimate the instantaneous phases of all the oscillators after the transformation of the attractors, which 
now yields monotonically increasing phase. Open circles in Figs. [IJb-f ) indicate the Poincare section. 

In the absence of coupling (e = 0.0) between the piecewise linear time-delay systems in the array, the 
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Fig. 5. (a) Time averaged phase (0^) given by ([7} and (b) time averaged frequency (n,) given by ([8]) of all the systems in the 
array Q with open end boundary conditions plotted as a function of the system index i for different values of the coupling 
strength e. 

individual subsystems evolve independently in an asynchronous fashion which is indeed confirmed from the 
monotonous increase in the phase differences, Atpi = (pi — (f>i, between the drive and the uniformly selected 
response systems (i = 10, 20, 30, 50) indicated by asterisk symbols in Figs. El Phase slips (marked by filled 
squares and triangles) in Figs. [3] for intermediate coupling strengths indicate that the coupled systems are 
in their transition to phase synchronization with the drive. For suitable values of e, the response systems 
in the array are entrained to the drive one sequentially, indicated by filled circles, as seen in Figs. El To 
be more precise, it is clear from this figure that, the 10th oscillator is synchronized first at e = 0.45, while 
the other systems away from it are still in their transition to phase synchronization. Further increase in e 
to e = 0.55 results in synchronization of the 20th oscillator with the drive leaving the remaining far away 
oscillators from it in their transition state. The 30th and 50th oscillators are entrained to the drive at 
e = 0.62 and 0.79, respectively, illustrating the existence of GPS through the sequential synchronization 
of the subsystems locally in the array. 

The existence of GPS is also confirmed by the average phase difference, r]{t), defined as 



v{t) 



N 



1 ^ 



;)• 



(6) 
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Fig. 6. Oscillator index vs oscillator index snap shots indicating the sequential phase synchronization and the dynamical 
organization of cluster states for different values of coupling strength. The different shapes (colors) indicate that the corre- 
sponding nodes are phase synchronized, (a) non-phase synchronized case for e = 0.0 (b) First fourteen oscillators in Eq. (jlJj) 
are phase synchronized with the drive system for e — 0.45. (c), (d) and (e) Sequential phase synchronization and the formation 
of small cluster states for e — 0.51, 0.58 and 0.62, respectively, and (f) global phase synchronization for e = 0.8. 



A global measure to characterize the existence of GPS is the average phase difference (f?(t)), which is 
shown in Fig. |4] for different values of e. The monotonous increase in r/(t) for e = 0.0, correspond to the 
independent evolution of all the subsystems in the array. Phases of the nearby oscillators are mutually 
locked sequentially as e is gradually increased, which is indicated by the low degrees of ri{t) for e = 0.56 
and 0.77 in Fig. |4l Finally, r/(t) ~ at e = 0.79 corroborating the existence of GPS in the array. 

The nature of dynamical organization of local oscillators to form GPS can be better understood by 
examining the time average phase and the time average frequency as a function of the oscillator index for 
different values of the coupling strength. For this purpose, we define the time averaged phase {{(pi)) as 



{(j)i{t)) = {2TTk + 2TT 

and the time averaged frequency as ((f^j)) 



t' - 1\ 

fi _ - 



lim — 

T^ooT 



(4 < i* < t\+i)i 



<j)i{t)dt, 



(7) 



(8) 



where t^ is the time of the k^^ crossing of the flow with the Poincare section of the i^^ attractor and {■..)t 
denotes time average. 

Figures [5)^a) and EJ^b) depict the average phase and the average frequency, respectively, for different 
e as a function of the oscillator index i. The random distribution of ((/)j(t)) (Fig. [Ha)(i)) and (^Qi{t)^ 
(Fig. [Hb)(i)) for e = 0.0 corresponds to the large difference in the initial frequency of the uncoupled 
oscillators due to the nonlinear parameter mismatch. The scenario of sequential phase synchronization to 
reach GPS can also be clearly visualized using the snap shots of the oscillators in the index vs index plot 
as shown in Fig. [H The diagonal line in Fig. ^a) for e = 0.0 clearly shows the asynchronous evolution of 
all the oscillators in the array. 

The phase locking of the first 14 oscillators with the drive, forming the main cluster, for e = 0.45 is 
shown in Figs. [Ha)(ii) and [5jb)(ii), while some of the other oscillators (oscillators with indices 22 — 24, 
35 — 36 and 41 — 45) away from it form phase synchronized clusters among themselves. The main cluster 
and other small clusters of phase synchronized oscillators are much more clearly visualized in Fig. [D^b) 
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for the same coupling strength, in which oscillators with same frequency are assigned the same shape and 
color. The average phase and the average frequency for e = 0.58 depicted in Figs.[5]^a)(iii) and Fig. [5]||b)(iii) 
indicate that the first 26 oscillators are synchronized with the drive with a few other small synchronized 
clusters away from it. Similarly, the first 35 oscillators are synchronized with the drive for e = 0.62 (see 
Figs. [5][|a)(iv) and Fig. [Hb)(iv)) along with other small clusters. Thus, it is clear that an increase in 
the coupling strength results in desynchronization of some of the existing clusters and formation of new 
clusters away from the main cluster. At the same time, nearby desynchronized oscillators from the main 
cluster becomes synchronized sequentially with it thereby increasing its size, contributing to the basic 
mechani sm of format i on of GPS via sequential synchronization of local oscillators as explained in detail 
in Ref. Suresh et all |2010(]. Finally for e = 0.8, all the oscillators become phase locked to attain GPS 
as depicted in Figs. [5l^a)(v) and [5|^b)(v), which remains stable for further larger values of the coupling 
strength. For the above values of e, the main cluster and other small clusters are also clearly seen in the 
index vs index plot in Figs. [6l 

Further, we have also plotted the average phase {{(f>i)) and the average frequency {{0,i)) of all the N 
oscillators as a function of e in Figs. [7| to get a global picture of the above dynamics. Broad distribution of 
{(j)i) and (Oj) for small e correspond to asynchronous states while for e > 0.75 both the average phase and 
average frequency converge to a single value indicating the existence of GPS. However, for the intermediate 
values of e, the main cluster and other small clusters are hardly visible in Figs. [71 Nevertheless, these clusters 
are clearly seen in Figs. [5] and [U where the snapshots of the dynamical organization of all the N oscillators 
in the array are depicted for specific values of the coupling strength. 

The frequency difference (Aili^j = Qi — ilj) and the frequency ratio (Qi/Qi) for the oscillators with 
the index i = 10, 20, 30, 40 and 50 along with their average (represented by black filled circles) of all the 
(A^ — 1) response oscillators as a function of the e are depicted in Fig. [8^a) and Fig. \Slh), respectively. 
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Fig. 8. (a) The frequency diflterence (A^i^^, i = 10, 20, 30, 40 and 50) and (b) the frequency ratio (f^i/^i, i = 10, 20, 30, 40 
and 50) of uniformly selected systems are plotted as a function of the coupling strength e G (0, 1.0). Each line corresponds to 
the difference/ratio between the response system and the drive system. The black filled circles indicate the average frequency 
difference/ratio of all the (A'' — 1) response systems from the drive system. 



It is also clearly seen from these figures that the nearby oscillators are synchronized sequentially as e is 
increased, that is, i = 10 is synchronized first and then i = 20, and so on as e is increased. Finally, all the 
N = 50 oscillators in the array are phase synchronized for e > 0.75 attributing to the emergence of GPS 
via local sequential synchronization. 



The well-known Kuramoto order parameter [Moreno &: Pachecd . l2004l | can also be used for quantifying 
the global phase synchronization in the array, which is defined as 



Re 



tip 



1 ^ 

n1^ 



j(t>j(t) 



(9) 



Here 4'j{t) represents the instantaneous phase of the j^^ system, tpit) is the average phase and (...)( denotes 
a time average. When all the systems are in a phase synchronized state, the value of i? ~ 1.0. Figure [9] 
shows that R decreases at first for small increase in e as observed in Figs. El^a) andElJb) accounting for the 
increased randomness in the phases of the coupled systems as seen in the average phase and the average 
frequency of all the oscillators in Figs. [71 Further increase in e results in a gradual increase in the value of 
R, reaching i? ~ 1.0 for e > 0.75 confirming the existence of GPS in the array of coupled piecewise linear 
time-delay systems. 

It is also to be noted that the phenomenon remains qualitatively the same even when the total number 
of oscillators in the array is increased. We also wish to emphasize that the results remain qualitatively 
unaltered even for different sets of random values for the nonlinear parameters, a^, confirming the robustness 
of our results. 
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Fig. 9. Plot of the phase order parameter (R) as a function of the coupling strength indicating global phase synchronization 
in the array of (A'^ — 50) coupled piecewise hnear time-delay systems Q. 



3.3. Confirmation of GPS from the original untransformed 
non-phase-coherent attractor 



Next, we use t wo d ifferent approaches, namely recurr ence quantification analysis [Romano et all l2005l : 
Marwan et all l2007l | and the concept of localized sets [Pereira et all l2007l | without estimating the phase 
explicitly to prove the existence of GPS from the original nonphase coherent hyperchaotic attractors. 



3.3.1. Recurrence quantification for GPS 

Now, we use the recurrence quantification measures discussed in Sec. 12.21 in characterizing the existence 
of GPS via sequential synchronization. Figure [10] depicts the generalized autocorrelation function of the 
drive Pi{t) and that of some of the response systems Pio{t), Psoit), and P^oii) for different e. For e = 0.0, 
none of the maxima of Pi{t) coincide in Fig. [TOTa) indicating asynchronous state. Pi{t) for i = 1, 10,30 
and 50 are plotted in Figs. fTOT b) and [TOlfc) for e = 0.45 and e = 0.8, respectively. In Fig. [TOlfb). the 
maxima of Pi{t) and Pio{t) are in perfect agreement (see Fig. llOlfblfi)! confirming the existence of phase 
synchronization between the oscillators with the index i = 1 and i = 10, which is in agreement with the 
results of the earlier section. On the other hand, some of the maxima of Pi(t) and Pso^t) in Fig.[TOlfb)(ii) are 
in agreement attributing to the transition to GPS among them and none of the maxima of Pi{t) and P^o^t) 
in Fig. [Torb)(ii) are in agreement corresponding to independent evolution, which confirms the sequential 
synchronization. All the maxima of Pi{t) are in complete agreement in Fig. fTOTc) confirming the existence 
of GPS. Also, the formation of clusters by the other asynchronous oscillators in the array can be confirmed 
by plotting their respective generalized autocorrelation functions, which will show that all their maxima 
are in good agreement with each other, whereas there exists a drift betwe en them and the maxima of the 
sequentially synchronized cluster as discussed in Ref . [Suresh et all l2010l | . 

The global scenario of the existence of GPS via sequential phase synchronization can also confirmed 
using the index GPR in analogy with the nature of the Kuramoto order parameter R, the average frequency 
difference and the average frequency ratio. The index GPR of the response systems i = 10, 20, 30, 50 
with that of the drive is shown in Fig. [TT] as a function of e. It is clear that the nearby oscillators to 
the drive are synchronized first as the coupling strength is gradually increased contributing to the local 
sequential synchronization. The results are exactly similar to that observed from the average frequency 
difference and the average frequency ratio in Figs. [8] indicating the phenomenon of GPS via sequential 
phase synchronization. 
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Fig. 10. Plots of the Generalized autocorrelation functions of the drive Pi{t) and some selected response systems (i = 
10, 30, and, 50) Pio{t), -P3o(i), and P^oit) indicating (a) Non-phase-synchronization for e = 0.0, (b) generahzed autocorrelation 
functions for e = 0.45 (bi) PS between the systems 1 and 10, (bii) approximate PS between the systems 1 and 30 and (biii) 
non PS between the systems 1 and 50, and (c) PS between all the systems (i — 1, 10, 30, and 50) for e = 0.8. 
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Fig. 11. Plots of the index CPR as a function of the coupling strength e. Different lines correspond to the CPR of different 
(i — 10, 20, 30 and 50) response systems with the drive system. The filled circles correspond to the mean value of the CPR of 
all the {N — 1) piecewise linear systems in the array. 



3.3.2. GPS using the concept of localized sets 

Next, we use the framework of localized sets described in Sec. 12.31 to demonstrate the existence of GPS 
via sequential phase synchronization. We defined the event as a Poincare section chosen at x{t — r) = 0.8 
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Fig. 12. First row (a-q) corresponds to the attractors of the drive system (i — 1) and the rows (b-r, c-s, and d-t) correspond 
to the attractors of some selected response systems (i = 10,30,50). In (a-d) the sets (represented by the filled circles) are 
spread over the attractors and hence there is no CPS for the value of coupling strength e — 0.0. In (e-h) for e — 0.4 and in 
(i-1, m-p, q-t) the sets are localized confirming the existence of GPS in the array for e — 0.52, 0.65, and 0.8, respectively. 



and x{t) < 0.8 on the attractor. The sets (indicated as filled circles) obtained by observing the drive 
system (i = 1) whenever the defined event occurs in the response system i = 10 along with the drive 
attractor are depicted in the first row in Fig. [T2j Similarly, the sets obtained by observing the response 
systems (i = 10, 30 and 50) whenever the event occurs in the drive system are depicted along with the 
corresponding response attractors in the other subsequent rows. 

The sets spread over the entire attractors in Figs. [T2lfa-d) for e = 0.0 confirm the asynchronous 
evolution of the subsystems in the array. The sets localized in a large area of the attractors in Figs. [T2re-f) 
for e = 0.4 indicate the transition to CPS among the corresponding oscillators, whereas the other oscillators 
away from it still evolve independently as indicated by the spread of the sets over their attractors (see 
Figs. [T2r g-h)). The more localized sets in Figs. [T2lfi-j) for e = 0.52 indicate that the respective oscillators 
are in complete phase synchrony while the other systems away from it are in asynchronization with it as 
confirmed by the spread of the sets in Figs. [T2rk-1) . For e = 0.65, the sets are bounded to smaller regions 
of the attractors in Figs. [EJ^m-o) attributing to the existence of CPS between the drive and subsystems 
10 and 30. The other subsystem with index i = 50 is in its transition state as indicated by the localized 
set but in a large region of the attractor as seen in Fig. [T2lfp). Thus it is clear from these figures that the 
oscillators away from the drive system are synchronized sequentially as the coupling strength is increased 
gradually. Finally, for e = 0.8, the sets are localized to a narrow region on the attractors (Figs. [I2jq-t)) of 
all the subsystems thereby confirming the existence of GPS. 
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Fig. 13. The first six maximal Lyapunov exponents Xmax of the threshold piecewise linear time-delay system (Eq. ((TD}) for 
the parameter values a — 1.2, /3 — 1.0, r G (1, 10) in the absence of the coupling e. 

4. Global Phase Synchronization in an Array of Piecewise Linear Systems with 
Threshold Nonlinear Function 

In this section, we demonstrate the existence of the GPS in another piecewise hnear time-delay system 
with a threshold nonlinear function. This system has been studied recently for its hyperchaotic nature 
even for small values of time-delay [see for details, Laks hmanan &: Senthilkum ar, 2010] and has also been 
experimentally realized using analog electronic circuits |Srinivasan et all |2011|] . Very recently, the chaotic 
phase synchronization has been experimentally confi rmed in this piecewise linear time-delay system along 
with numerical simulation Senthilkumar et all l201Cll | and various types of synchronization t r ansiti ons have 



been demonstrated both experimentally and numerically in this system [Srinivasan et all 1201 ll ]. In this 
paper, we consider a linear array of piecewise linear systems as in Eq. ^ and the nonlinear function f{x) 
is now chosen to be a piecewise linear function with a threshold nonlinearity, 



f{x) = AF* - Bx. 



(10) 



Here 



p* 



— X 
X, 



-X < x'^ 

-X < X < X* 
X > X 



(11) 



where x* is a controllable threshold value and A and B are positive parameters. This function f{x) employs 
a threshold controller for flexibility. It effectively implements a piecewise linear function. The control of 
this piecewise linear function facilitates controlling the shape of the attractors. Even for a small delay 
value, this system exhibits hyperchaos and can produce multi-scroll chaotic attractors by just introducing 
more threshold values. In our analysis, we chose x* = 0.7, A = 5.2, B = 3.5, ai = 1.2, /? = 1.0, r = 6.0 
and the nonlinear parameters Oj, of the response systems in the array are chosen randomly in the range 
at G (1-18, 1.24). Note that for this set of parameter values a single uncoupled system exhibits a highly 
complicated hyperchaotic attractor with three positive Lyapunov exponents. The first six largest Lyapunov 
exponents of the uncoupled system are shown in Fig. [13] as a function of time-delay r E (1, 10). 

We have not yet succeeded in generalizing the nonlinear transformation (Eq. ([I])) to capture the phase 
of the non-phase-coherent hyperchaotic attractor of the above systems due to the multiscroll attractor. 
However, we find that the recurrence-based indices serve as excellent quantifiers in identifying the transition 
from non-synchronized to phase synchronized state both quantitatively and qualitatively. We have also 
characterized the occurrence of GPS using the concept of localized sets. 
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Fig. 14. Generalized autocorrelation functions of the drive Pi{t) and some of the selected response systems {i = 10, 30, and, 
50) Pio{t), P3o(i), and P5o(i) indicating (a) Non-phase-synchronization for e = 0.0, (b) generalized autocorrelation functions 
for e — 0.8 (bi) PS between the systems 1 and 10, (bii) approximate PS between the systems 1 and 30 and (biii) non-PS 
between the systems 1 and 50, and (c) PS between all the selected systems (i — 1, 10, 30, and 50) for e — 1.2. 



4.1. GPS using recurrence analysis 

We have calculated the CPR and the generalized autocorrelation function P{t) to confirm the existence 
of the GPS in the array. The generalized autocorrelation function of the drive Pi{t) and that of some 
response systems {i = 10,30, and 50), Pio{t), Psoit), and P5o{t), are depicted in Fig. [lUfor different values 
of the coupling strength. In the absence of coupling (e = 0.0), all the systems evolve independently and 
hence the maxima of their respective generalized autocorrelation functions do not occur simultaneously 
as shown in Fig. [T^a). On increasing the coupling strength, the oscillators with a lower value of index 
in the array become synchronized first resulting in sequential phase synchronization and this can also be 
identified from the generalized autocorrelation functions of the response systems in the array. For instance, 
Pio(t), P3o(i), and P5o{t) are shown along with Pi{t) in Fig. [HT b) for e = 0.8. It is clear from this figure 
that the maxima of the drive Pi{t) and those of the response Pio{t) are in complete agreement with each 
other [Fig. [Hlfb)(i)] indicating the existence of PS between them. On the other hand, only some of the 
maxima of the response system Psoit) are in coincidence with those of the drive [Fig. [T^ b)(ii)] illustrating 
that the response system i = 30 is in transition to PS, whereas the maxima of the response system P5o{t) 
do not coincide with those of the drive [Fig. [T^ b)(iii)] indicating that the response system i = 50 is in 
an asynchronous state for the same value of e. For e = 1.2, almost all of the positions of the peaks of the 
generalized autocorrelation functions Pi(t), Pio(t), Psoit), and P^oit) are in agreement with each other as 
illustrated in Fig. [HTc) confirming the existence of GPS via sequential phase synchronization. It is also to 
be noted that the magnitudes of the peaks of all the oscillators have generally different values and these 
differences in the heights of the peaks indicate that there is no correlation in the amplitudes of the coupled 
systems. 

The dynamical organization of GPS via sequential synchronization and the clustering can be visualized 
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Fig. 15. Snap shots of the node vs node diagrams indicating the sequential phase synchronization and the organization 
of cluster states for different values of coupling strength for system @ with the threshold nonlinearity (|lUp . (a) non-phase 
synchronized case for e = 0.0 (b) First four oscillators are phase synchronized with the drive system for e — 0.65. (c), (d) and 
(e) Sequential phase synchronization and the formation of small cluster states for e = 0.7, 0.94 and 1.13, respectively, and (f) 
global phase synchronization for e = 1.25. 

clearly by plotting some snap shots of the oscillators in the index vs index plots as shown in Fig. [T5j The 
oscillators that evolve with identical phase are assigned with identical shapes (colors). The diagonal line 
in Fig. [T5lfa) for e = 0.0 correspond to the oscillator index i = j and the oscillators evolve independently. 
Figure [T5lfb) indicates that the first four oscillators in the array are synchronized with the drive for 
e = 0.65 along with the five small separate clusters. In Fig. [TSTc) for e = 0.7 the first six oscillators form 
a synchronized cluster along with eight other small clusters. Similar small clusters are shown in Fig. [T5lfd) 
and (e) for e = 0.94 and e = 1.13, respectively, in addition to the single large cluster formed by sequential 
phase synchronization. Finally GPS of all the systems in the array is illustrated in Fig. [TSTf) for e = 1.25. 
The existence of GPS via sequential phase synchronization is also quantified using the index CPR ^ 
of the response systems with the drive as shown in Fig. [TBI The different lines correspond to the index 
of the oscillators {i = 10, 20, 30, 50) in the array. It is evident from the figure that the oscillators with 
increasing index attain the value of unity in a sequence as a function of the coupling strength and finally 
for e > 1.15 the CPR of all the response systems with the drive reaches unity confirming that all the 
coupled oscillators are in GPS. The mean value of CPR of all the response systems in the array, shown as 
filled circles, also confirms the existence of GPS for e > 1.15. 



4.2. 



GPS using the concept of localized sets 



We have also confirmed the existence of GPS in the linear array of threshold piecewise linear time-delay 
systems (Eq. (fTO]) ) by using this concept of localized sets as in the previous case. Now, we will demonstrate 
the existence of GPS via sequential phase synchronization in some randomly selected response systems {i 
= 1, 10, 30, 50). The set obtained by sampling the time series of one of the systems whenever a maximum 
occurs in the other system is plotted along with the attractors of the same systems. The set, indicated as 
filled circles, obtained by observing the drive system (i = 1) whenever the maxima occurs in the response 
system (i = 10) is shown in Fig. fTTT a) and that obtained by observing the response systems i = 10, 30, 50 
whenever the maxima occurs in the drive system are shown in Figs. [TTT b-d) for the value of coupling 
strength e = 0.0. As the sets are spread over the attractors, all the systems evolve independently and there 
is no CPS in the absence of coupling between them. Further when we increase the coupling strength to e 
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Fig. 16. The index CPR as a function of the coupling strength e. Different lines correspond to the CPR of different (i = 
10, 20, 30, and 50) response systems with the drive system. The filled circles correspond to the mean value of the CPR of all 
the {N — 1) piecewise linear systems in the array. 

= 0.7, the oscillator {i = 10) is partially synchronized with the drive as the sets are almost localized but 
the sets in the oscillators i = 30 and 50 are spread over the attractor which means that they are not yet 
phase synchronized with the drive system. This is shown in Figs. [T71fe-h). Again increasing the coupling 
strength to e = 0.85, the sets are further bounded to a small region over the attractors which shows that 
the oscillator i = 10 is synchronized with the drive, but the oscillator i = 30 is partially synchronized where 
the sets are almost localized and i = 50 is not yet phase synchronized with the drive as represented by the 
spread of the sets over the attractor in Figs. [TTT i-l). Further, the Figs. [T71fm-p) and Figs. fTTT q-t) indicate 
the situation for e = 1.0 and e = 1.2, respectively, where all the oscillators are now phase synchronized 
with the drive as the sets are localized over the attractor confirming the existence of GPS in an array via 
sequential phase synchronization as the coupling strength is increased. It is to be noted that the sets are 
not yet completely localized in Figs. [TTjj for e = 1.0, whereas it is localized for e = 1.2 (Figs. 1171). 



5. Partial Phase Synchronization (PPS) in a Linear Array of Piecewise Linear 
Time-delay Systems with Ring Topology 

So far we considered the dynamics of coupled arrays of time-delay systems with open end boundary 
conditions. Now, we wish to demonstrate the existence of a partial phase synchronization (PPS) in an array 
of time-delay systems with a ring topology. In this coupling configuration, all the systems are not fully 
synchronized, instead they split into several subgroups to form different phase synchronized clusters. The 
phenomenon may also be called cluster synchronization. Here the systems within ea ch cluster rnaintain 



perfect phase synchr onization. This k ind of PPS mostly occ urs in neural net works Denzl et alx . 120081 ] 



coupled chaotic systems [Bjorn et al. 



chemical oscillations [Kiss et all 12001 1 and El-Nino systems [Stein et all l201l[ |. and has been studied in 



20061 1 as well. Also the transition from nonsy nchronization to pha, se 



Zhuang et all . l2002l |. 



synchronization via PPS has been studied in two-dimensional coupled map lattices 

The dynamical equation of a linear array of coupled time-delay systems with closed end boundary 
conditions is given as 



Xi{t) = -I3xi{t) + aif{xi{t - r)) + e(xj_i(t) - Xi{t)), 



(12) 



where i = 1,2,- ■ ■ , N, with the following periodic boundary conditions: xq = xn and xn+i = xn- The 
function f{x), and the system parameters have been chosen as in Sec. 13.11 In this type of coupling config- 
uration the signal of the A^th system is fed into the first system and so there is no specific drive system 
and that each system sends its signal to the nearby system in a unidirectional way. In the absence of the 
coupling, there is no synchronization among the systems. If we increase the coupling, the systems with the 
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Fig. 17. First row (a-q) corresponds to the attractors of the drive system (i = 1) and the rows (b-r, c-s, and d-t) correspond 
to the attractors of some selected response systems (i = 10,30,50). In (a-d) the sets (represented by the filled circles) are 
spread over the attractors and hence there is no CPS for the value of coupling strength e = 0.0. In (e-h) for e = 0.7 and in 
(i-1, m-p, q-t) the sets are localized confirming the existence of GPS in the array for e = 0.85, 1.0, and 1.2, respectively. 



same phase/frequency form a small group of clusters. If we increase the coupling above a threshold value 
{£thr)i PPS occurs with a relatively large group of phase clusters. For even larger values, one finds that the 
PPS state collapses. 

The nonlinear transformation used in Sec. [3] is no longer valid to transform the non-phase-coherent 
attractor into a phase-coherent attractor due to the new boundary condition. However, we have already 
found that the recurrence-based indices serve as excellent quantifiers in identifying phase synchronization 
in non-phase-coherent attractors both qualitatively and quantitatively. Further, we also characterized the 
occurrence of the cluster formation using the concept of localized sets. We will use these quantifications in 
the following. 

We have calculated the CPR of each of the {N — 1) systems in the array by taking z = 1 as the 
reference system (but we also confirmed that the dynamics does not change even if we take any other 
system in the array as a reference system). In Fig. [18] we have plotted the CPR of every system in the 
array as a function of the system index (z) for different values of coupling strength. In the absence of the 
coupling (e = 0.0), the systems in (fT2]l evolve independently as indicated by the random distribution of 
CPR in Fig. [TSTa). Increasing the coupling strength to e = 2.22 results in synchronous evolution of a few 
of the oscillators, leading to groups of small clusters (which are having the same CPR values) as seen in 
Fig. [TST b). However there is no single global synchronized state in the present network in contrast to the 
sequentially synchronized single cluster as discussed in Sec. El Above a threshold value {ethr > 3.0), we can 
observe the occurrence of a PPS in the array with large groups of phase synchronized cluster. Figs. flST c-f) 
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Fig. 18. The index CPR as a function of the system index i for various values of couphng strength e in ring topology for 
(fT2)) . (a) For e = 0.0, (b) for e = 2.22, (c) for e = 3.08, (d) for e = 3.14, (e) for e = 3.72 and (f) for e = 4.0. 
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Fig. 19. Snap shots of the node vs node plots for various values of coupling strengths showing the occurrence of phase clusters, 
(a) For e = 0.0, (b) for e = 2.22, (c) for e = 3.08, (d) for e = 3.14, (e) for e = 3.72 and (f) for e = 4.0. 



indicate the occurrence of the PPS with separate groups of phase synchronized clusters for e = 3.08, 3.14, 
3.72 and 4.0, respectively. If we increase the coupling strength for even larger values, PPS continues to 
exist until e increases to very large values (e > 50) where it collapses. For very large values of the coupling 
strength the clusters lose their stability resulting in a desynchronized array. 

The formation of the phase synchronized clusters can also be clearly illustrated and visualized by 
using the snap shots of the oscillators in the index vs index plot as shown in Fig. [T9l In this figure the 
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Fig. 20. The index CPR as a function of the couphng strength e. Different hnes correspond to the CPR of different (i — 
10, 20, 30, 40 and 50) systems. The fiUed circles correspond to the mean value of the CPR of all the (A*' — 1) piecewise linear 
systems in the array. 
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Fig. 21. The occurrence of phase clusters is also explained using the concept of localized sets, (a) For e = 0.0 figure shows 
the attractor of the system i = 30 and the sets of the systems i = 17, 18, 19,22 and 31, (randomly chosen) (b) for e — 3.72 
with the sets of the systems inside the cluster i = 17, 28, 31, 34 and 37 (which lie within the cluster), (c) for e — 3.72 with the 
sets of the systems outside the cluster i — 18, 19, 22 and 27. 



systems with identical phase (having the difference in their CPR value less than 0.001) are assigned with 
identical shapes (colors). The diagonal line in Fig. [TWa) for e = 0.0 correspond to the oscillator index 
i = j and the oscillators evolve independently. Fig. [TWb) displays the occurrence of small groups of phase 
synchronized clusters in distant nodes (i = 2 and 50) for e = 2.22. If the coupling is increased further the 
nodes with the same frequency form groups of clusters with each cluster having different phases resulting 
in PPS. In contrast, in Fig.[6]due to the open end boundary conditions, the nearest subsystems to the drive 
start to synchronize with it along with the formation of small groups of clusters. Further increase in the 
coupling strength results in the decomposition of the small clusters which then synchronize with the main 
cluster. Consequently GPS results in the array with open end boundary conditions. Figs. [TWc-f) show the 
occurrence of PPS in the network with the formation of large phase clusters for e = 3.08,3.14, 3.72 and 
4.0 respectively. 

The mechanism for the formation of partial phase synchronization can be explained as follows: In 
the absence of the coupling the individual systems oscillate independently with different phases due to 
the mismatches in the nonlinear parameters a^. Upon increasing the coupling strength, due to the ring 
configuration, every individual oscillator starts to drive the nearest oscillator and the oscillators with small 
differences in their phases/frequencies in the array synchronize themselves to form small groups of clusters 
leaving the other oscillators with large phase differences to evolve independently. Further increase in the 
coupling strength results in an increase in the sizes of the clusters due to the locking of the phases of the 
nearby oscillators to the cluster for appropriate e and finally ending up with a large cluster near the center 
of the ring for large e. However there exists other small clusters in the ring contributing to PPS and hence 
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Fig. 22. The index CPR as a function of the system index i for various values of coupling strength e in the array (|12p with 
ring topology of the threshold piecewise linear system, (a) For e = 0.0, (b) for e = 6.12, (c) for e = 6.80, (d) for e = 10.4, (e) 
for e = 11.05, and (f) for e = 11.8. 
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Fig. 23. Snap shots of the node vs node plots for various values of coupling strengths shows the occurrence of the phase 
clusters in threshold piecewise linear systems, (a) For e = 0.0, (b) for e = 6.12, (c) for e = 6.80, (d) for e = 10.4, (e) for 
£ = 11.05, and (f) for e = 11.8. 

there is no GPS in such an array. As we have seen earher in Sec. El in GPS, as long as the couphng increases, 
the nearest subsystems of the drive start to synchronize with it and the remaining asynchronous systems 
with the closest frequencies form small groups of phase clusters. Further increase in the coupling results in 
the formation of a single large cluster with the drive due to the decomposition of the small clusters with 
different frequencies and thereby resulting in GPS. But in the ring configuration, if the coupling increases 
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Fig. 24. The index CPR as a function of the couphng strength e of the threshold piecewise hnear time-delay system. Different 
lines correspond to the CPR of different {i — 10, 20, 30, 40 and 50) systems. The filled circles correspond to the mean value of 
the CPR of all the {N — 1) piecewise linear systems in the array. 

further (above the threshold value) the distant nodes with same phase/frequency form groups of clusters 
with each cluster having its own phase resulting in PPS. 

The synchronization transition to PPS is quantified using the value of CPR as a function of the 
coupling strength. In Fig. [20] the different lines correspond to the CPR of the randomly selected systems 
in the network (i = 10, 20, 30, 40 and 50). The filled circles correspond to the mean value of the CPR of all 
the (A^ — 1) piecewise linear systems. It is evident from this figure that only very few systems get phase 
synchronized to form small clusters and globally there is no phase synchronization when the coupling is 
below the threshold value (sthr < 3.0) which is indicated by the mean value of CPR near to zero. Once 
the coupling strength crosses the threshold value (e ~ 3.0), there is a sudden jump to CPR ~ 0.8 which 
is an indication of the occurrence of PPS. It is evident from this figure that the system i = 50 is already 
synchronized with the drive reference system xi{t) as it is one of the nearest neighbors, which is also clearly 
seen in Figs. [18] and [T9j One may note the interesting fact that this transition is of first order (sharp) type, 
while in the case of open end boundary condition (see Sec. [3]) it is smooth and of second order type. 
The reason for this kind of transition is as follows: In open end boundary conditions, the subsystems are 
sequentially synchronized with the drive, so the GPS increases gradually in the array as a function of the 
coupling strength. Hence the transition is continuous (second order). But in a ring topology all the systems 
attain PPS suddenly at a particular value of the coupling strength due to the mutual sharing of the signals 
which leads to the discontinuous (first order) transition. 

Next we calculate the localized sets by defining the event among the systems within the cluster and 
observing the other systems that are in the same cluster. In this case the obtained sets of the systems within 
the cluster are localized on the attractor, while the sets of the systems outside the cluster gets spread over 
the entire attractor. For illustration, in Fig. [2]T a) the attractor of the system i = 30 is plotted along with 
the sets of the randomly selected systems i = 17, 18, 19, 22 and 31 (which are represented by different 
shapes and colors) where the sets are spread over the attractor which corresponds to the unsynchronized 
state for e = 0.0. Fig. [2lT b) corresponds to e = 3.72 for the sets of the systems 17,28,31, and 37 which 
are all within the cluster. The figure shows that the sets of the corresponding systems are localized on the 
attractor in the same place implying that these systems form a clustered state. On the other hand, for 
the same value of the coupling, we have plotted the sets of the systems outside the cluster (i = 18, 19, 22 
and 27). This shows that the sets are spread over the entire attractor indicating that the corresponding 
systems are not in phase synchronization (Fig. [2T]^c)). 



6. Partial Phase Synchronization in an Array of Piecewise Linear Systems with 
Threshold Nonlinearity 

In this section, we demonstrate the existence of the PPS in the second piecewise linear time-delay system 
with a threshold nonlinear function (jlOp . The form of the piecewise function /(x) in Eq. (J12p and the 
system parameters are set to be the same as in Sec. [5] 
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Fig. 25. The occurrence of the phase clusters in threshold piecewise linear systems in the ring topology is explained using the 
concept of localized sets, (a) For e = 0.0 shows the attractor of the system i = 20 and the sets of the systems i = 10, 15, 19, 18, 24 
and 40, (b) for e = 11.05 with the sets of the systems inside the cluster i = 15, 16, 18 and 26. (c) for e = 11.05 with the sets 
of the systems outside the cluster i = 10, 13, 24 and 40. 



In Fig. [22j the CPR is shown as a function of the system index (i) for different values of the couphng 
strength. For e = 0.0, in Fig. [22]^a), there is no correlation in the values of CPR of the systems attributing 
to the desynchronization state. Figures [22]^b) andl22lfc) correspond to the coupling strength e = 6.12, and 
6.80, respectively, which display the occurrence of small phase synchronized clusters. If we increase the 
coupling strength e > 10.0 one can observe the occurrence of PPS in the array with large groups of phase 
synchronized clusters. Figures [22lfd-f) indicate the occurrence of the PPS with separate groups of phase 
synchronized clusters for e = 10.4, 11.05 and 11.8, respectively. 

The dynamical organization of the cluster formation is clearly visualized in Fig. 1231 node vs node plot, 
for various values of the coupling strength. The different shapes (colors) correspond to the nodes which 
are in the different phase synchronized clusters. In Fig. [23]^a) the diagonal line shows the oscillator index 
i = j evolving independently for e = 0.0. Figures [23lfb) and [23]^c) display the occurrence of several small 
phase synchronized clusters for e = 6.12 and 6.80, respectively. If the coupling is increased further, the 
nodes with the same frequency form groups of clusters with each cluster having different phases resulting 
in PPS as shown in Fig. [23lfd-f) for the values of e = 10.4, 11.05 and 11.8, respectively. 

The transition from nonsynchronization to partial phase synchronization is quantified by plotting the 
CPR as a function of the coupling strength. In Fig. [23] the different lines correspond to CPR of the randomly 
selected systems (i = 10, 20, 30, 40 and 50) and the filled circles represent the mean value of the CPR of 
(N — 1) systems. It is evident from this figure that there is no synchronization (CPR< 0.2) below e < 10.0. 
There is a sudden jump (a first order transition) in the value of CPR (to ~ 0.85) for e ~ 10.0 indicating 
the onset of PPS in the array along with the formation of a large group of phase synchronized clusters. 

The occurrence of these phase synchronized clusters in the array of threshold piecewise linear time- 
delay systems is also demonstrated using the concept of localized sets. In Fig. \25\ we have plotted the 
attractor of the system i = 20 along with the sets of some randomly selected systems in the array. Each 
shapes (colors) correspond to the sets of different subsystems. Figure. [25ja) is plotted for e = 0.0 and the 
sets of the systems i = 10, 15, 18, 24 and 40 spread over the entire attractor indicating the asynchronized 
state. For e = 11.05, the sets of the systems (i = 15, 16, 18 and 26) inside the large cluster are localized on 
the attractor as shown in Fig. [25r b) illustrating that the respective systems are in a phase synchronized 
state. On the other hand, for the same value of the coupling strength, we have plotted the sets of the 
systems outside the cluster {i = 10,13,24 and 40), which are spread over the attractor i = 20 indicating 
that these subsystems are not in a phase synchronized state (see Fig. 125(c)). 



7. Summary and Conclusion 

We have demonstrated the existence of global and partial phase synchronizations in an array of unidirec- 
tionally coupled nonidentical piecewise linear time-delay systems with two different boundary conditions. 
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Coupled with our earlier work on Mackey-Glass system Suresh et aUl201Cll ]. our studies clearly establish 



the generic nature of the underlying phenomena. In particular, in a linear array with open end boundary 
conditions, we have demonstrated the emergence of GPS via sequential synchronization. We have shown 
that in addition to the main phase synchronized cluster centered at the drive, the remaining asynchronous 
systems from the main cluster organize themselves to form different clusters for low values of the cou- 
pling strength. Further increase in the coupling strength leads to the formation of a single large cluster 
resulting in GPS by a decomposition of the other clusters. The synchronization transition is of second 
order type. We have confirmed the existence of GPS via sequential phase synchronization by estimating 
the phase difference, the average frequency and the average phase as a function of the oscillator index and 
the coupling strength from the transformed attractors. Furthermore, we have also confirmed the existence 
of GPS from the original non-phase-coherent attractors of the coupled piecewise linear time-delay systems 
using appropriate recurrence quantification measures and the concept of localized sets without explicitly 
estimating the phase. We have also confirmed the existence of GPS via sequential phase synchronization 
in a second piecewise linear system with threshold nonlinear function. 

The existence of a partial phase synchronization (PPS) is demonstrated in an array with closed end 
boundary conditions (ring topology) and the synchronization transition is of first order type. PPS is 
corroborated using recurrence analysis and the concept of localized sets using the original non-phase- 
coherent hyperchaotic attractor of the piecewise linear time-delay systems. Further, the mechanism for the 
formation of both GPS and PPS are elaborated. We have also confirmed the existence of PPS in an array of 
Mackey-Glass time-delay systems with closed end boundary conditions; however the results are not given 
here as they are similar in nature to the above systems. 

This study can also be extended to two-dimensional lattices of coupled time-delay systems with and 
without delay coupling and in networks of time-delay systems with indirect global dynamic environment 
coupling, as well as networks with complex topology. 
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